The idea of using Josephson junctions as sources of 
electromagnetic radiation is promising owing to their 
small dimensions, good tunability, and capability of op- 
erating at frequencies up to several hundred gigahertz 0| . 
However, the power of radiation available from a single 
junction is not sufficient for most applications, which ne- 
cessitates using arrays of junctions. A stack arrangement 
of long Josephson junctions (LJJs) ||^ arouses interest 
because of possible improvements of the properties of 
LJJ oscillators in terms of impedance matching, output 
power, and integration level. The low- Tc superconductor 
thin film technology allows growth of high-quality mul- 
tilayers with many Josephson tunnel barriers (for exam- 
ple, Nb/Al-AlOx/Nb) stacks). Moreover, the discovery 
of an intrinsic Josephson effect in some high- Tc super- 
conductors such as BSCCO convincingly showed that 
these materials are essentially natural superlattices of 
LJJs formed on the atomic scale Recent theoretical 
investigations and experiments showed that the induc- 
tive coupling between adjacent junctions leads to diverse 
and nontrivial dynamic behaviour patterns of such struc- 
turesglll. 

One possible way to produce coherent radiation from 
stacks of LJJs is to form a regular Josephson vortex lat- 
tice (JVL) and move it by an external current. In order 
to produce maximum radiation from a stack at a given 
lattice velocity and external magnetic field, a rectangu- 
lar arrangement of vortices, when vortices in neighboring 
layers are located one over another, is most preferable. 
Such a lattice is feasible provided the corresponding so- 
lution is stable. In the present paper we show that in 
a two-junction stack a rectangular JVL is unstable at 
low velocities, and stability of such a solution can be 
achieved at high velocities of JVL, provided there is a 
large external magnetic field, or large enough damping, 
or small stack length. Specifically, the above conditions 
may explain the results of numerical ^ ^ and exper- 
imental investigations in which a possibility of exis- 
tence of rectangular JVL in stacks of two and more LJJs 
is shown. Our result complements the results reported 
in 1^ where the authors argue that a rectangular JVL 
is stable at high velocities in an infinite stack. The in- 
completeness of this result is explained by the fact that 
the authors of |^ have not taken into account all possible 
perturbations of the solution in the form of a rectangular 
JVL. 

Let us consider the set of equations describing a sim- 
plest layered structure — a stack consisting of two LJJs 
with magnetic coupling between the layers ||, |[ ^: 
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[d"^ + jdt)ipi^2 + sin ipi^2 - j- (1) 



Here 931,2, 7, j are the Josephson phase difference, damp- 
ing constant, and bias current density, respectively. The 
magnetic coupling parameter is denoted by 77. It is de- 
termined by the formula 

7? = A(dsinh-|- + 2Acoth-^)-\ (2) 



where A is the London penetration depth, t is the thick- 
ness of the superconducting layer, d is the distance be- 
tween two superconducting layers. We start with the 
assumption that the system is infinite in space. 

To investigate a two-junction stack it is convenient to 
introduce new variables ip± — (ipi ± ip2)/'2, which obey 
the set of equations: 

cj.dlip± = {dt + ^dt)Lp± + sin Lp± cos ^p^ - j±, (3) 
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where = 1, — {1 ~ rj) / {1 + rf), j+ = 
have renormalized the coordinate Xnew 
in Eqs. d). 

The set of Eqs. (||) has a solution describing rectan- 
gular JVL. Assuming the external magnetic field to be 
high we can write down the analytical expressions for 
this solution [foi: 



(^1 = h{x — ut) + Im- 
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where L — —h?{l — u^) -f- i^uh, u is the JVL veloc- 
ity, /i 3> 1 is the dimensionless external magnetic field. 
Velocity u and damping 7 are related through the en- 
ergy balance condition [ pH which is actually the current- 
voltage characteristic of the stack with the rectangular 
JVL: ^ ^ 

-^uh ~\ — Im— . (5) 
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We note that the solution in the form is valid only 
provided \L^^ \ <C 1. If ^17 1 the previous condition is 
satisfied at all velocities otherwise it breaks in a region 
near u = 1. This region corresponds to the peak in the 
current- voltage characteristic (^). 

In order to investigate the stability of rectangular JVL 
we search for the solution of (||) in the form Lp± = 
f^. -t- 5'p± where tf \ are given by (Q) and 5ip± are small 
perturbations (|5</'±| <C 1). Substituting this solution 
into (H) and neglecting the terms nonlinear in 5lp± we 
obtain: 
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i^^J^ w cos/i(x- ut) +Re((l - e2'''("^-"*))/2L). 
We will refer to 5lp+ and Sip^ as symmetrical and anti- 
symmetrical perturbations, respectively. The set (|^) is 
actually two independent equations so we can analyze 
them separately. Thus we have divided the problem of 
rectangular JVL stability into two ones — the problems 
of stability with respect to symmetrical and antisymmet- 
rical perturbations. 

We start with an analysis of the "subluminal" (u < 1) 
solution stability. The problem of a rectangular JVL sta- 
bility with respect to symmetrical perturbations is sim- 
ilar to the problem of stability of the periodical vortex 
chain in LJJ which was solved in 0. Thus a "sublu- 
minal" rectangular JVL is stable to symmetrical pertur- 
bations. To analyze the stability of with respect to 
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antisymmetrical perturbations we use Eq. (^) for Sip^. 
This equation is a relativistic invariant with c_ being the 
characteristic velocity of antisymmetrical perturbation. 
From (H) and the expression for c_ it is seen that c_ < 1 
at any stack parameters. In other words, antisymmet- 
rical perturbations are always slow compared with the 
symmetrical ones. Therefore, to investigate the "sub- 
luminal" solution stability it is necessary to distinguish 
between two cases: u < c_ and u > c_. 

Let us first consider the case u < c_. We perform the 
Lorentz transformation in Eq. (||) for Stp-: 



vt 



t — {v/c^_)x 



with velocity v — u. Introducing ip = 5ip^ we obtain the 
equation 
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-.{ipr — Ulj^l) + COS ip]^ ■ -0 



where the parameter depends only on the coor- 
dinate ^. After the renormalization of the co- 
ordinate h^Jl — {u/ c^)'^^oid — "i^new and time 
h^Jl — {u/ C-Ytou = 2t„p,„ and introduction of the 
small parameter fi = 



4/i ^(c?_ — u^) ^ we have: 
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ij^^+T^tp^ = — (V'TT+^rV'T)+M[cos2^-|-Re- 
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where — F,- — fih'yc^/2. We look for 

the solution of in the form of the Fourier integral 

= r 

Fourier image of iP{Ctt) is 



ip{^,uj)e-''^^^. The equation for the 
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According to the Bloch theorem, the solutions of 
have the form tp{^) = exp (ig^)J7g(^), where q is the 
quasimomentum and Uq{^) is the function with the pe- 
riod TT. Let us find the eigenfrequency spectrum Lo{q). 
At /i = the spectrum is 



(9) 



At /X ^ the maximum perturbations of the spectrum 
are achieved near the middle point (q — 0) and the edges 
(q = ±1) of the first Brillouin zone. In the vicinity of 
g = we search for the solution of (S) in the form 
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where ao,a±2 are constants. After substituting it 
into (^) we obtain the dispersion characteristic as the 
condition for ao,±2 at which the solution of (^) is not 
equal to zero: 



(10) 




0,0 0,2 0,4 



1 ,0 1 ,2 



Figure 1: Eigenfrequency uj{q) spectrum as a function of 
quasimomentum q for antisymmetric perturbations (case 

M < c_) 

where a = Re(2L)~^ — 7^/4. Near q = 1 the solution 
has the form 



V;(e) =e*('?"i)«[aie'«-|-a_ie^ 



where a±i are constants. Substitution of this expression 
into (ph gives 
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Far from the middlepoint and the edges of the first Bril- 
louin zone the spectrum remains unperturbed and is 
given by the formula (||). 

The dependencies of real and imaginary parts of eigen- 
frequency uj{q) are shown in Fig.|^. It is seen that at 
small q some roots of the dispersion equation ( p^ have 
positive imaginary parts. This means that perturbations 
with small q will exponentially grow with time, i.e. the 
solution (^) is unstable. We see that the instability which 
is obvious in the case of low density chains (due to repul- 
sion of vortices in the neighbouring layers) is not changed 
by stability in the case of denser chains. This result is in 
agreement with the one obtained in and reflects the 
fact that the vortex chains in the neighbouring layers 
tend to shift and form the triangular JVL. 

Let us now consider the case of high velocities c_ < 
M < 1. As before, we perform the Lorentz transformation 
in Eq. (^ for Sip- but now with the velocity v = c^_/u. 
Introducing tjj = 6ip^ we obtain the equation 
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Mr ^V'5)+COS(^J -V", 



where the parameter depends only on the time r. This 
equation turns into 



V'rT + TrV'T = c?_(V'j5 +r^V^)-A^(cos2T-|-Re 
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(12) 
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Figure 3: Stability diagram for the rectangular vortex 
lattice (case hj < 1). Latin letters mark the regions of 
different behaviour of small perturbations: a - longwave 
instability, b - shortwave instability to antisymmetric 
perturbations, c - shortwave instability to both sym- 
metric and antisymmetric perturbations, d - stability. 
Hatching shows the region where the solution (Q) is not 
valid. 



Figure 2: Quasienergy s{k) as a function of wavenumber 
k for antisymmetric perturbations (case c_ < u < 1) 
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fi = Ah'^{u^ - cl)-\ = /iM/17/2, Fj = /x/17/2. We 
look for the solution of ( |l2|) in the form of the Fourier 
integral ■0(^, t) = ip{k, T)e~*''^ ||. The equation for 
the Fourier image of r) is 
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-iP^^+Tri^r = -c^_{k'^ -ikT^)xp- ^j,{cos2T+Re^-—- — 

(13) 

The solution of ( p^ ) has the form iI^{t) = 
exp (— ier)C/£(r) where e is the quasienergy and 
Us{t) is the function with the period tt. Let us find the 
spectrum e(fc) by the scheme used above. In the vicinity 
of £ = the spectrum is 
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and near e = 1 
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It follows from the Eq. (|T^) that at certain k the roots 
of the dispersion equation may have positive imaginary 
part as it is shown in Figj^. Ime(fc) reaches its maximum 
value at fc = cl^. At small 7 {lTii£{k))max ~ — Fr/2. 
Therefore, the perturbation with k 
time like 



depends on 



^{k,T) ~ exp[(^ - 
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It is seen from this expression that the solution (^ is 
either stable or unstable depending on sign of the differ- 
ence — u-f. If huj < 1 the solution (||) is paramet- 
rically unstable. The region of k corresponding to the 
growing perturbations is equal to 



Afc ; 
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This parametric instability may be "suppressed" either 
by increasing the external magnetic field h or by increas- 
ing the damping 7. We would like to emphasize that the 
instability appears due to the periodicity of the solution. 
Therefore the results obtained in for the isolated 
vortices cannot be applied to the case of periodic vortex 
chains. 

It remains to investigate the stability of the "super- 
luminal" JVL {u > 1). We start with the analysis of 
the stability with respect to symmetrical perturbations. 
Substituting the solution (^) into Eq. (H) for S(p+ and 
applying the Lorentz transformation 
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vr 



with V = u ^, we obtain 
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V'gC=V'rr-<- , =(-0r--V'c)+COSy^O ■V', (17) 

^l-(l/u)^ U 

where ^ = S^p^ . In this equation the coefficient at V' de- 
pends only on time r and is periodical. Analyzing this 
equation by the method used for the case C- < u < 1 
we arrive at the conclusion that there arc shortwave per- 
turbations which depend exponentially on time t and 
the solution (^) may be stable or unstable depending on 
sign of the difference h^^ — wy. The same result is ob- 
tained also for the antisymmetric perturbations in the 
case M > 1. The main feature of this result is that the 
stability condition h^^ < uj may be achieved by increase 
of u. However, the alternate component of the electric 
field ipt of such a solution has the order h~^u{u'^ — 1)~^ 
and is small. Therefore, the power of radiation from the 
edge of the stack in this case is not sufficient for appli- 
cations in spite of the fact that the JVL is rectangular. 

Combining the results obtained above we build the 
diagram (Fig. ^) showing the regions of stability and 
instability in terms of the lattice velocity u. At < u < 
c_ the solution (^ is unstable with respect to longwave 
perturbations and formation of the triangular JVL. At 
c_ < ii < (hj)^^ the solution is unstable due to the 
parametric resonance and the instability growth rate is 
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proportional to —wy. In the region near u — 1 where 
the condition \L^^\ ^ 1 breaks the solution (||) is not 
valid and more thorough investigation is required. Then, 
at M > {hj)~^ the parametric resonance is "suppressed" 
and the solution becomes stable. The diagram Fig. ^ is 
built for the case /17 < 1. When ^,7 > 1 the instability 
region is < u < {h'f)^^. At last, when {hj)~^ < c_ 
the instability region is < m < c_ . 

Consider now the case when the system is not infinite 
in space. The simplest way to change to a finite system 
is to set periodic boundary conditions 

(pi,2ix = l,t)^ ViAx = 0, <) + 27riV, 

where N is the number of vortices trapped in each junc- 
tion of the stack, / is the length of the system. The 
boundary conditions for perturbations are written as be- 
low: 

5V±{ineio = hl/2,T) = 5ip±{£,new = 0, t). (18) 

At w < c_ the conditions (|l^) lead to discreteness of 
quasimomentum q with the step Air/hl. As the insta- 
bility interval is located near q = 0, mode with q = 
will always be in this interval, i.e. it will grow with 
time. Hence, the periodic boundary conditions cannot 
provide stability of a rectangular JVL at small veloci- 
ties. At c_ < u < 1 the wavenumber k is discrete with 
the step An/hi. At large I when Afc > in/hi at least 
one of the possible values of k will be in the interval of 
instability (p^), therefore, the solution will be unstable 
at > u-f. The necessary condition of stability is 

I > nhc^{c^_ - u^). 

This means that to provide a rectangular JVL stability 
at c_ < u < 1 it is necessary to decrease the system 
length I so that no permitted value of k is in the in- 
terval (p^. At 1 < u < {h^)^^ the situation is more 
complicated because the instability intervals in k have 
different widths and are differently located for symmet- 
rical and antisymmetrical perturbations. But it is clear 
that the decrease of / lowers the number of permitted k 
in the instability intervals and may lead to stability. 

Let us summarize the obtained results. At low veloci- 
ties a rectangular JVL in an infinite two-junction stack is 
unstable with respect to perturbations with scales much 
larger than the lattice period. The physical meaning 
of this is that the vortex chains in the first and sec- 
ond junctions tend to shift with respect to each other to 
compose a triangular JVL. This is in agreement with the 
result of . At high but " subluminal" velocities the lat- 
tice may be unstable with respect to perturbations with 
scales comparable to the period of the JVL. The growth 
rate of this parametric instability is determined by the 
difference —uj and, therefore, the instability may be 
"suppressed" by high enough external magnetic field, or 
by the damping in the system, or by the increase of the 



lattice velocity. The existence of this instability com- 
plements the results obtained in j|] where the authors 
did not take into account the short-wave perturbations. 
In the present paper we also show that the stability of 
rectangular JVL at low damping is possible in a finite 
two-junction stack. Thus we argue that the formation of 
a rectangular JVL which is reported in @, ||, |[ is asso- 
ciated either with the finiteness of the system or with the 
suppression of the instability by three ways mentioned 
above. 
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